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Abstract. We review our contribution to infrared Renormalization Group (RG) effects to General 
Relativity in the context of galaxies. Considering the effective action approach to Quantum Field 
Theory in curved background, we argued that the proper RG energy scale, in the weak field hmit, 
should be related to the Newtonian potential. In the galaxy context, even without dark matter, this 
led to a remarkably small gravitational coupling G variation (about or less than 10^'^ of its value 
per hght-year), while also capable of generating galaxy rotation curves about as good as the best 
phenomenological dark matter profiles (considering both the rotation curve shape and the expected 
mass-to-light ratios). Here we also comment on related developments, open issues and perspectives. 
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1. INTRODUCTION 

Currently there is a large body of data coming from cosmological and astrophysical ob- 
servations that is mostly consistent with the existence of dark matter. Such observations 
also suggest that the hypothesized particles that constitute dark matter have very small 
cross section and travel at speeds much lower than Ught. These lead to the cold dark mat- 
ter framework, which is one of the pillars of the current standard cosmological model 
ACDM. 

It is not only tempting, but mandatory to check if such dark matter particles exist 
(by detecting them in laboratory based experiments, for instance) and also to check if 
the gravitational effects that lead to the dark matter hypothesis could follow from a 
more detailed and complete approach to gravity. The effects of pure classical General 
Relativity at galaxies have been studied for a long time and, considering galaxy rotation 
curves, the differences between General Relativity and Newtonian gravity are negligible, 
since in a galaxy matter moves at speeds much lower than that of light and is typically 
subject to weak gravitational fields (4> <^ c^), which leads to the Newtonian limit of 
General Relativity \ 



There are some proposals that consider General Relativity in the context of galaxies which do not lead 
to Newtonian gravity, see e.g. [1, 2, 3, 4, 5]. It is not impossible that a reasonable explanation for galaxy 



There is however a newer approach that may change considerably the role of dark 
matter, while following standard physical principles. Namely, the investigation of the 
running of the gravitational coupling parameter G on large scales as induced by the 
renormalization group framework. 

The running of coupling constants is a well known phenomenon within Quantum 
Field Theory. It is well-known that the renormalization group method can be extended to 
quantum field theory on curved space time and to some models of quantum gravity (see, 
e.g., [6]), such that the beta functions can be interpreted in this framework. Concerning 
the high energy (UV) behavior, there is hope that the running of G in quantum gravity 
may converge to a non-Gaussian fixed point (asymptotic safety) [7, 8]. Our present 
concern is, however, not about the UV completeness, but with the behavior of G in 
the far infrared regime (IR). In the electromagnetic case the IR behavior corresponds 
to the Appelquist-Carazzone decoupling [9] (see e.g., [10] for a recent derivation of 
this theorem). In the case of gravity the same effect of decoupling has been obtained in 
[11, 12], but only for the higher derivative terms in the gravitational action. It remains 
unclear whether such decoupling takes place for the other terms. This possibility was 
studied on phenomenological grounds a number of times before, e.g. [13, 14]. 

In [15] we presented new results on the application of renormalization group (RG) 
corrections to General Relativity (GR) in the astrophysical domain. Previous attempts 
to apply this picture to galaxies have considered for simplicity point-like galaxies. 
We extended previous considerations by identifying the proper renormalization group 
energy scale ji and by evaluating the consequences considering the observational data 
of disk galaxies. Also we propose a natural choice for the identification of jU, linking 
it to the local value of the Newtonian gravitational potential. With this choice, the 
renormalization group-based approach is capable to mimic dark matter effects with 
great precision. This picture induces a very small variation on the gravitational coupling 
parameter G, namely a variation of about 10~^ of its value across 10^ light -years. We call 
our model RGGR, in reference to renormalization group effects in General Relativity. 

In order to evaluate the observational consequences of the RGGR model and to 
compare it to other proposals, recent high quality observational data [16, 17] from 
nine regular spiral galaxies were mass-modelled using the standard procedures for the 
baryonic part, and four different models for the "dark" component: i) the RGGR model; 
ii) one of the most phenomenological successful dark matter profiles, the Isothermal 
profile [18]; and two alternative models which were built to avoid the need for the dark 
matter: iii) the Modified Newtonian Dynamics (MOND) [19, 20] and iv) the Scalar- 
Tensor- Vector Gravity (STVG) [21]. The latter is a recent relativistic proposal that is 
capable of dealing with galaxy rotation curves and other phenomena usually attributed 
to dark matter. For galaxy rotation curves phenomena, STVG becomes equivalent to a 
similar proposal called MSTG [22, 23, 24]. The model parameters that we use to fit 
galaxies in the STVG framework can be found in ref. [23]. While for MOND we use the 
ao value as given in [25]. 

The quality of the rotation curve fits and total stellar mass as inferred from the RGGR 



rotation curves may rely on similar approaches, nevertheless up to now none of such proposals have found 
a baryonic mass distribution that is in conformity with astrophysical expectations. 



model is perfectly satisfactory considering both the general behavior of the model and 
its results when applied to nine particular galaxies, as analyzed in [15]. It is about the 
same of the Isothermal profile quality, while it seems significantly better than the quality 
of MOND and STVG. In the case of MOND, we did numerical experiences with ao as 
a free parameter and found that, albeit the concordance with the shape of the rotation 
curve can considerably increase in this case, the concordance with the expected stellar 
mass-to-light ratios remains unsatisfactory (similar conclusions have also appeared in 
some recent papers, e.g. [17, 26, 27], and it seems that the concordance can only be 
improved by adjusting the MOND's ju(;c) function in an ad-hoc way). 



2. THE RUNNING OF G 

The j8 -function for the gravitational coupling parameter G has been discussed in the 
framework of different approaches to Quantum Gravity and Quantum Field Theory in 
curved space-time. In [15] we followed the derivation used previously in [13]. If G does 
not behave as a constant in the far IR limit, it was argued in [13] (and recently in more 
details in [28]) that the logarithmic running of G is a direct consequence of covariance 
and must hold in all loop orders. As far as direct derivation of the physical running of G 
is not available, it is worthwhile to explore the possibility of a logarithmically running 
G at the phenomenological level. 

Consider the following infrared j8 -function for General Relativity, 

Equation (1) leads to the logarithmically varying G(ju) function, 

where jUq is a reference scale introduced such that G(/io) = Gq. The constant Go is the 
gravitational constant as measured in the Solar System (actually, there is no need to be 
very precise on where G assiraies the value of Gq, due to the smallness of the variation 
of G). The dimensionless constant V is a phenomenological parameter which depends 
on the details of the quantum theory leading to eq. (2). Since we have no means to 
compute the latter from first principles, its value should be fixed from observations. It 
will be shown that even a very small v can lead to observational consequences at galactic 
scales. 

The action for this model is simply the Einstein-Hilbert one in which G appears inside 
the integral, namely. 




In the above, G should be understood as an external scalar field that satisfies (2). Since 
for the problem of galaxy rotation curves the cosmological constant effects are negligi- 
ble, we have not written the A term above. Nevertheless, for a complete cosmological 
picture, A is necessary and it also runs covariantly with the RG flow of G (see e.g., [13]). 



There is a simple procedure to map solutions from the Einstein equations with the 
gravitational constant Go into RGGR solutions. One need not to follow this route, one 
may find all the dynamics from the RGGR equations of motion, which can be found by a 
direct variation of the action (3) in respect to the metric, leading to equations of motion 
that have the same form of those of a scalar-tensor gravity^. In this review, we will 
proceed to find RGGR solutions via a conformal transformation of the Einstein-Hilbert 
action, and to this end first we write 

G = Go + 5G, (4) 

and we assume 5G/Gq <^ 1, which will be justified latter. Introducing the conformally 
related metric 

81J.V = -^8nv, (5) 

the RGGR action can be written as 

S,Ug]=SM+0{5G^). (6) 

where S^,, is the Einstein-Hilbert action with Gq as the gravitational constant. The above 
suggest that the RGGR solutions can be generated from the Einstein equations solutions 
via the conformal transformation (5). Indeed, within a good approximation, one can 
check that this relation persists when comparing the RGGR equations of motion to the 
Einstein equations even in the presence of matter [15]. 

In the context of rotation curves of galaxies, standard General Relativity gives essen- 
tially the same predictions of Newtonian gravity. The Newtonian potential 4>n.„, is related 
to the metric by 

^00 = -(1 + "^). (7) 

Hence, using eq. (5), the effective RGGR potential <I> in the non-relativistic limit is given 
by 

c2 5G 

^ = ^.^+-—. (8) 
2 Go 

An equivalent result can also be found from the evaluation of a test particle geodesies 
[15]. In the context of weak gravitational fields ^>n„,/c^ <^ 1 (with 4>n,„ = at spatial 
infinity) holds, and hence the term 5G/Gq should not be neglected. 

In order to derive a test particle acceleration, we have to specify the proper en- 
ergy scale jU for the problem setting in question, which is a time-independent gravi- 
tational phenomena in the weak field limit. This is a recent area of exploration of the 
renormalization group application, where the usual procedures for high energy scat- 
tering of particles cannot be applied straightforwardly. Previously to [15] the selection 
of jU 1/r, where r is the distance from a massive point, was repeatedly used, e.g. 
[29, 30, 31, 32, 13]. This identification adds a constant velocity proportional to v to 



We stress that it is only the from since RGGR is not a type of scalar-tensor gravity, and G is not a 
fimdamental field of the model. 



any rotation curve. Although it was pointed as an advantage due to the generation of 
"flat rotation curves" for galaxies, it introduced difficulties with the TuUy-Fisher law, 
the Newtonian limit, and the behavior of the galaxy rotation curve close to the galactic 
center, since there the behavior is closer to the expected one without dark matter. In [15] 
we introduced a /i identification that seems better justified both from the theoretical and 
observational points of view. The characteristic weak-field gravitational energy does not 
comes from the geometric scaling 1/r, but from the Newtonian potential 4>n„„,. However, 
the straight relation ji <^ 4>n„, leads to /i 1 /r in the large r Umit; which is unsatisfactory 
on observational grounds (bad Newtonian limit and correspondence to the TuUy-Fisher 
law). One way to recover the Newtonian limit is to impose a suitable cut-off, but this 
does not solves the TuUy-Fisher issues [13]. Another one is to use [15] 



where ^>o and a are constants. Apart from the condition < ^>o < (i.e., essentially 

4>o is a reference Newtonian potential), the precise value of $o is largely irrelevant for 
the problem of rotation curves. The relevant parameter is a. It is a phenomenological 
parameter that depends on the mass of the system, and it must go to zero when the 
mass of the system goes to zero. This is necessary to have a good Newtonian Umit. 
From the TuUy-Fisher law, it is expected to increase monotonically with the increase 
of the mass. Such behavior is indeed found from the galaxy fits done in [15]. In a 
recent paper, an upper bound on vet in the Solar System was derived [28]. In galaxy 
systems, va|oj„, ~ 10~^, while for the Solar System, whose mass is about 10-10 of that 
of a galaxy, va|s„,„,sy„„„ < 10^^^. It shows that a linear decrease on a with the mass is 
sufficient to satisfy both the current upper bound from the Solar System and the results 
from galaxies. 

We also point that the above energy scale setting (9) was recently re-obtained from a 

more theoretical perspective [33]. 

Once the ji identification is set, it is straightforward to find the rotation velocity for a 
static gravitational system sustained by its centripetal acceleration. 



Contrary to Newtonian gravity, the value of the Newtonian potential at a given point 

does play a significant role in this approach. This sounds odd from the perspective of 
Newtonian gravity, but this is not so from the General Relativity viewpoint, since the 
latter has no free zero point of energy. In particular, the Schwarzschild solution is not 
invariant under a constant shift of the potential. 

In the following, we will comment on the effect of the relation (10) to galaxy rotation 
curves. First from a more general perspective, and then to the modeling of individual 
galaxies. 




(9) 




(10) 



3. GALAXY ROTATION CURVES 



Before proceeding to specific galaxy rotation curves modeling, it is more instructive to 
analyze general features of the relation (10), and to compare it to the standard approach. 
In [15] we analyze some general aspects and scaling laws of the RGGR model, with 
no dark matter, in comparison to the isothermal profile; both of them, at this step, 
without gas and with an exponential stellar disk. In particular, it was pointed that the 
RGGR rotation curves have a reasonable shape to fit galaxies (i.e., no clear problems 
like increasing or decreasing too fast, oscillations...), and that they effectively behave 
similarly to cored dark matter profiles at inner radii, whose effective core radius scales 
with the galaxy disk scale length. Further details in our paper. 

We have also extended the previous analysis by adding a gas-like contribution (a re- 
scaled version of the NGC 3198 gaseous part). In particular, this numerically evaluates 
how the model behaves on the presence of density perturbations at large radii. In the first 
plot of fig. (1) it is displayed the result for RGGR, which is remarkably good (a similar 
plot can be found in our original paper), while in the others plots in fig. (1) (presented at 
the Conference, but not in [15]) one sees the results for the same mass distribution but 
different choices for the energy scale^ jl. 




FIGURE 1. The additional circular velocity squared induced by different choices of jx. The first 
plot refers to RGGR, the two others to different identifications of jj.: one depends on the Newtonian 
acceleration (a variation inspired on MOND) and the other on the (baryonic) matter density. All the plots 
above display the additional squared velocity of each model divided by = vac^ and as a function of 
R/Rd, where R is the the radial cylindrical coordinate and Rd is the stellar disk scale length. Black lines 
depict the additional velocity due to a pure exponential stellar disk, while the gray solid lines take into 
account the gas mass Mgas for different values of / = Mgas/M^^^^, with / = 0.2,0.7, 1.2, ...,9.7 (i.e., 
the black lines stand for / — 0). See [15] for further details. 

From fig. (1), the two other proposals different from RGGR are seen to be unsuited 
as replacements for dark matter. In particular, both are too sensitive to the gas presence, 
and both eventually add negative contributions to the total circular velocity at large radii. 

In [15] we used a sample of nine high quality and regular rotation curves of disk 
galaxies from [16, 17]. In figs. (2, 3) we show one of ours results (see [15] for the 
complete set and further details) in comparison to the results of three other models: 
a cored dark matter profile (Isothermal profile), the Modified Newtonian Dynamics 
(MOND) and the recently proposed Scalar-Tensor- Vector Gravity (STVG). 



^ These other choices are also imsatisfactory from the theoretical perspective, since they have no direct 
relation to the local energy of the gravitational field in the weak field regime (^Newt ^ c^)- 




FIGURE 2. NGC 2403 rotation curve fits. The red dots and its error bars are the rotation curve 
observational data, the gray ones close to the abscissa are the residues of the fit. The solid black line 
for each model is its best fit rotation curve, the dashed yellow curves are the stellar rotation curves from 
the bulge and disk components, the dotted purple curve is the gas rotation curve, and the dot-dashed green 
curve is the resulting Newtonian, with no dark matter, rotation curve. 




FIGURE 3. Stellar disk mass-to-light ratio in the 3.6jUm band as a function of the color J — K. 
Each galactic disk is represented above by an open circle, with a reference error bar of 50% of the Y^,o 
value. The black open squares display the 7*0 values and their associated la errors for each galaxy as 
inferred from the rotation curve fits for each model. The highlighted square and circle correspond to the 
NGC 2403 galactic disk mass-to-light ratios. See [15] for further details. 



Due to the considerably large uncertainty in the total stellar mass of each stellar 
component (disk and bulge), we first use the total stellar mass as a free parameter for the 
fittings (achieved from a minimization considering the errors). At a second stage, we 
compare the resulting value with stellar population expectations, following the standard 
approach. 



On the free parameters of each model, we remark that besides the total stellar mass, 
the Isothermal profile has two additional free parameters, the RGGR model has a single 
free parameter (a) while MOND and STVG have no free parameters that can vary from 
galaxy to galaxy. On the other hand, both of the latter depend on constants whose values 
are calibrated considering its best fit in a large sample of galaxies. We remark that the v 
parameter in RGGR cannot vary from galaxy to galaxy, but a can, and galaxy rotation 
curves are sensible to the combination va, whose value is about the order of 10~^. The 
best fit for NGC 2403 yields va = (1.66± 0.01) x 10-"^. 

4. CONCLUSIONS 

We presented a model, motivated by renormalization group corrections to the Einstein- 
Hilbert action, that introduces small inhomogeneities in the gravitational coupling across 
a galaxy (of about 1 part in 10^) and can generate galaxy rotation curves in agreement 
with the observational data, without the introduction of dark matter as a new kind of 
matter. Both High and Low Surface Brightness galaxies were tested [15] . Considering 
the samples of galaxies evaluated in [15], the quality of the RGGR rotation curves, to- 
gether with the corresponding mass-to-light ratios, is about the same than the Isothermal 
profile quality, but with one less free parameter. We expect that similar results would 
hold in regard to other cored dark matter profiles, while our results seem better than 
those achieved by the NFW profile [15]. We also compared the results of our model 
with MOND and STVG, and at face value our model yielded clearly better results. 

Our results can be seen as a next step compared to the previous models motivated 
by renormalization group effects in gravity, e.g. [13, 34]. Their original analyses could 
only yield a rough estimate on the galaxy rotation curves, since they were restricted to 
modeling a galaxy as a single point. Trying to extend this approach to real galaxies, 
we have shown that the proper scale for the renormalization group phenomenology is 
not of a geometric type, like the inverse of the distance, but is related to the Newtonian 
potential with null boundary condition at infinity. 

The essential feature for the RGGR rotation curves fits is the formula (10), which is 
by itself a simple formula that provides a very efficient description of galaxy rotation 
curves. 

There are several tests and implications of this model yet to be evaluated. In particular 
we are working on applying the RGGR framework to a larger sample of galaxies 
(including elliptical galaxies) [35] and galaxy-galaxy strong lensing [36]. Related work 
on CMB, BAO and LSS in search for a new cosmological concordance model is also a 
work in progress [37]. 
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